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ABSTRACT 



Context. The analysis of solar-like oscillations for stars belonging to a binary system provides a unique opportunity to probe the 
internal stellar structure and to test our knowledge of stellar physics. Such oscillations have been recently observed and characterized 
for the A component of the 70 Ophiuchi system. 

Aims. We determined the global parameters of 70 Ophiuchi AB using the new asteroseismic measurements now available for 70 Oph A 
and tested the input physics introduced in stellar evolution codes. 

Methods. Three different stellar evolution codes and two different calibration methods were used to perform a comprehensive analysis 
of the 70 Ophiuchi system. 

Results. A model of 70 Ophiuchi AB that correctly reproduces all observational constraints available for both stars is determined. An 
age of 6.2± 1.0 Gyr is found with an initial helium mass fraction Y\ = 0.266±0.015 and an initial metallicity (Z/X)[ = 0.0300±0.0025 
when atomic diffusion is included and a solar value of the mixing-length parameter assumed. A precise and independent determination 
of the value of the mixing-length parameter needed to model 70 Oph A requires accurate measurement of the mean small separation, 
which is not available yet. Current asteroseismic observations, however, suggest that the value of the mixing-length parameter of 
70 Oph A is lower or equal to the solar calibrated value. The effects of atomic diffusion and of the choice of the adopted solar mixture 
were also studied. We finally found that the different evolution codes and calibration methods we used led to perfectly coherent results. 

Key words, stars: oscillations - stars: interiors - stars: fundamental parameters - stars: binaries: visual - stars: individual: 70 Ophiuchi 



1. Introduction 

The solar five-minute oscillations have provided a wealth of in- 
formation on the internal structure of the Sun. These results 
stimulated various attempts to detect a similar signal on other 
solar-type stars by photometric or equivalent width measure- 
ments. In past years, the stabilized spectrographs developed 
for extra-solar planet search have finally achieved the accu- 
racy ne^de^Jo_deJe£^rolar ; like oscillations on other stars (see 
e.g. [Bedding & Kieldsedl2007h . A major difficulty is the con- 
frontation between these observations and theoretical models. 
The classical observational measurements available for an iso- 
lated star (effective temperature, metallicity, luminosity) com- 
bined with the oscillation frequencies provide strong constraints 
to the global parameters of the star but are often not sufficient to 
unambiguously determine a unique model and to really test the 
input physics included in the computation of the stellar mod- 
els (see for example the ca se of the isolated star f3 Virginis, 
Eggenberger & Carrier 2006). In this regard, binary stars consti- 
tute ideal asteroseismic targets in order to test our knowledge of 
stellar physics. In the case of a binary system, we can indeed as- 
sume that both stars share the same age and initial chemical com- 
position. The additional constraints imposed by this assumption 
are then extremely valuable to accurately determine the global 



properties of the stars. Moreover, in the case of a binary system, 
the masses of both components are accurately known by com- 
bining visual and spectroscopic orbits. 

Unfortunately, individual frequencies of solar-like oscil- 
lations have only been identified for stars belonging to 
two different binary syste ms: the G2 dwarf star a Cen 
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A jBouchy & Carried 12001 iBedding et aD l2004t iBlTzotetakl 
| 2007[) . the Kl dwarf star a Cen B dCarrier & Bourbanl 
\200A iKieldsen et al.l 120051) and the F5 IV-V star Procyon 
A dMartic et al.1 [2004 lEggenberger et alJ f2004a: M osser et all 
l2007h . Recently. ICarrier & E ggenbergenl2006l) reported the de- 
tection of solar-like oscillations on the bright K0 dwarf star 
70 Ophiuchi A (HD 165341A) belonging to the nearby visual 
binary system 70 Ophiuchi, which is among the first discovered 
binary systems. The aim of the present paper is to perform a 
comprehensive calibration of the 70 Ophiuchi system that in- 
cludes the new seismological data available for the A compo- 
nent. In addition to the determination of the 70 Ophiuchi AB 
global parameters, we also test and compare the theoretical tools 
used for the modeling of stars for which p-modes frequencies 
are detected by performing this analysis with three different stel- 
lar evolution codes and two different calibration methods. This 
work therefore takes place in the continuity of the work under- 
taken by the Evolution and Seismic Tools Activity (ESTA) group 
within the CoRoT mission (M onteiro et al. 2006). 



2 



P. Eggenberger et al.: Analysis of 70 Ophiuchi AB including seismic constraints 



Table 1. Orbital elements derived for the 70 Ophiuchi system 
when all available measurements (speckle and radial velocity) 
are included (see text for more details). 



Parameter 




Value 


P (days) 




32279 ± 15 


T (HJD) 




45809 ± 14 


e 




0.5005 ± 0.0006 


V (km s - 


) 


-7.026 ± 0.015 


i 




121.1 ±0.1 


id (deg) 




193.4 ± 0.3 


n (deg) 




121.7 + 0.2 


K i (km s" 


) 


3.51 ±0.04 


K 2 (km s~ 


) 


4.25 ± 0.05 


a (arcsec) 




4.526 ± 0.007 


Mi (Me) 




0.89 ± 0.02 


M 2 (Mo) 




0.73 ± 0.01 


it (mas) 




194.2 ± 1.2 



The observational constraints available for the 70 Ophiuchi 
system are summarized in Sect. 2. The input physics and the 
computational methods used for the calibrations are described in 
Sect. 3. The results are presented in Sect. 4, while the conclusion 
is given in Sect. 5. 

2. Observational constraints 

2.1. Binarity 

70 Ophiuchi is a very well studied stellar system composed 
of a K0 and K5 dwarf star with an orbital period of 88.38 yr 
dPour baix 2000). It is known as a binary by both spectroscopy 
and speckle interferome t ry. Ra dial- velocity curves were firs t 
obtained by iBatten et all d 19841) and IBatten & Fletcherl (Il99ll) . 
while the whole or b it (sp eckle and radial velocity) was pub- 
lished by iPour baix (2000). This binary system was observed 
over 6 nights in July 2004 with the Harps spectrograph mounted 
on the 3.6-m telescope at La Silla Observatory (ESO, Chili) 
dCarrier & Eggenberger! 120061) . Combining ou r new radial- 
veloc ity measurements with previous ones (Pourbai x et al.l 
2004), we use the ORBIT c ode made avai lable by T. Forveille 
and developed in Grenoble dForveille et al.ll 19991) to simultane- 
ously fit the radial velocities and the speckle measurements (see 
Fig.[TJand|2]i. Note that the very high precision of the new Harps 
data is quite important in order to derive an accurate value for K2 
and hence a reliable mass for the B component. Moreover, these 
new data add an important observation of the radial velocity of 
70 Oph B since only few measurements were previously ob- 
tained for this component (see Fig . [TJ . The values of the derived 
orbital parameters, of the parallax, and of the masses of both 
components are listed in Table [TJ We note that the value of the 
parallax derived in this way (194.2 +1.2 mas) is in g ood agree- 
ment w ith the value of 195.7+0.9 mas determined bv lSoderhielml 
( 1999) using the Hipparcos astrometry with ground-based obser- 
vations available at that time. 

2.2. Effective temperatures and chemical composition 

The effective temperature and iron abundance of the two stars 
were obtained from our own Harps spectra by f ollowi ng the 
methodology and line-lists used in Santo s etail d2004l) . This 
determination makes use of an iteractive procedure based on 




Phase 

Fig. 1. Radial velocity orbit of the 70 Ophiuchi system. Full 
and open symbols correspond to 70 Oph A and B, respectively. 
Triangles indicate the new Harps measurements, while circles 
indicate previous measurements. 

iron excitation and ionization equilibrium. We refer to this pa- 
per for more details. The results of this method were shown to 
give excellent results for solar-type stars, and the derived stel- 
lar parameters, namely the effective temperature, are compat- 
ible wit h those derived from o ther methods (e.g. the IR-flux 
method, Casa grande et al. I [20061) . For the metallicity of the sys- 
tem a mean value of [Fe/H] = 0.04 + 0.05 was thus determined 
and used for the present calibration. No a-element enhancement 
is suggested for the 70 Ophiuchi system by current spectroscopic 
data. Concerning the effective temperature, we derived a value 
of 5300 ± 50 K for the A component. For the B component, an 
effective temperature of 4205 K is found with an uncertainty es- 
timated to 100K. We however note that our spectroscopic de- 
termination is less adapted for this kind of cold star. We thus 
decided to adopt as T^ ff a mean value between our determined 
spectroscopic temperature of 4205 K, our photometric deduced 
value of 4220 K (see Sect ion 2.3) and the effect ive temperature 
of 4740 K determined by iLuck & Heiterl d2005i) . In this way a 
mean temperature of = 4390 K was determined for the B 
component. We also decided to adopt a large error of 200 K on 
the effective temperature of 70 Oph B, so that our calibration of 
the system is only weakly constrained by this value, which is not 
precisely determined. 

2.3. Luminosities 

From 1977 to 2005, the 70 Ophiuchi syst em has been mea- 
sured in the Geneva photometric system dGolavl [1 980) with 
the photoelectric photometer P7 dBurnet & Rufenerlll979l) in- 
stalled first on the Swiss telescope in La Silla (ESO, Chile) 
and then on the Belgian telescope at La Palma. Nine measure- 
ments were obtained. However, due to the small separation be- 
tween components and to too large a diaphragm, these pho- 
tometric measurements only represent the mean flux of both 
components V tot = 4.012 + 0.013, Vl tot = 4.779 + 0.013 and 
B2 tot = 5.323 + 0.015 mag. In order to obtain individual pho- 
tometric measurements of both components, the 70 Ophiuchi 
system was observed with the C CD camera C2 in stalled on 
the Swiss telescope at La Silla dCherix et al.1 12006) between 
2005 August 20-23 using very short exposure times of 0.02, 
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alpha [arcsec] 

Fig. 2. Apparent visual orbit of the 70 Ophiuchi system. The 
solid line indicates the line of nodes, while the square indicates 
the position of the periastron. 

Table 2. Observational constraints for 70 Oph A and B. 



70 Oph A 



70 Oph B 



n (mas) 




194.2 


+ 1.2 


M/M Q 


0.89 ± 


0.02 


0.73 ± 0.01 


V (mag) 


4.191 ± 


0.014 


6.061 ± 0.021 


L/L e 


0.53 ± 


0.02 


0.15 ±0.02 


7W (K) 


5300 


+ 50 


4390 ± 200 


[Fe/H] 




0.04 + 


0.05 


Av OuHz) 


161.8 


+ 0.8 





range 1 1 to 14 cm s" 1 . Unfortunately, the mono-site nature of the 
observations, coupled to the low resolution of the time series and 
to the faint signal-to-noise, do not allow a unique identification 
of the detected frequencies. Indeed, the small separation between 
1-2 and I — modes is not clearly detected and the £ = 2 modes 
cannot therefore be unambiguously disentangled from the i — 
modes. Consequently t = 0,2 can be replaced by i— 1 modes 
and {= 1 by 1-2 modes. Although these measurements can- 
not allow an unambiguous identification of individual detected 
oscillation frequencies, they provide us with a precise determi- 
nation of the mean large separation of 70 Oph A. By averaging 
the 4 different large sepa rations observed for the mode s listed as 
t = modes in Table 3 of lCarrier & Eggenberger! (120061) with an 
uncertainty on individual frequencies estimated to half the time 
resolution, we obtain a mean value Avo = 161.8+0.8 fiHz that is 
used in the present computations. We note that the error on this 
va lue is slightly larger than the uncertainty of 0.3 ^Hz quoted 
bv lCarrier & Eggenberger! d2006l) who included large separation 
averages taken between non-consecutive modes. 

Although the mean small separation of 70 Oph A cannot 
be clearly determined from current asteroseismic measurements, 
these observations nevertheless suggest that this value must be 
included in a frequency interval centered around the daily alias 
of 1 1 .57 fjHz introduced by the mono-site observations and de- 
limited by the frequency resolution of the time-series (2.2yuHz). 
Indeed, if the value of the small separation of 70 Oph A dif- 
fers from 11.57 ± 2.2 yuHz, we can then reasonably think that it 
should have been revealed by existing asteroseismic measure- 
ments. The inclusion of this additional observational constraint 
in the calibration of the 70 Ophiuchi system is discussed in detail 
in Sect. O 



3. Stellar models 



0.03 and 0.2 s in the V, VI and B2 bands, respectively. A 
large number of exposures were obtained in order to reduce the 
noise due in great part to the scintillation: 80, 18 and 12 ex- 
posures were thus taken in the V, VI and B2 filters, respec- 
tively. These CCD frames allowed us to determine the lumi- 
nosity shifts between both components; we thus found a dif- 
ference of 1.870 ±0.014, 1.920 ±0.014 and 2.255 ± 0.028 mag 
in the V, VI and B2 bands, respectively. Combining these 
values with the previous ones, we finally derived mean vi- 
sual magnitudes V A = 4. 191 ± 0.014 and V B = 6.061 ± 0.021 mag 
and photometric indexes (B2 - V1) A = 0.501 ± 0.020 and (B2 - 
V1) B = 0.836 ±0.040 for the 70 Ophiuchi system. Note that 
these indexes lead to an effective temperat ure of 4220 K for the 
B component (Ra mirez & Melendezll2005l) . 

Combining the mean magnitudes with the parallax deduced 
from the radial-velocity and speckle orbits {n- 194.2 ± 1.2), 
and by using the solar absolu te bolometric magnitude 
M bol>0 = 4.746 dLeieune et all Il998l) with the bolometric cor- 
rections from iFlowerl d 19961) (BC A = -0.20 ± 0.02 and BC B = - 
0.69 ± 0.17 mag), the luminosities are finaly determined with a 
value of 0.53 ± 0.02 L for 70 Oph A and 0.15 ± 0.02 L for 
70 Oph B. 

2.4. Asteroseismic constraints 

So lar-like oscillations in 70 Ophi uchi A were recently observed 
by ICarrier & Eggenberger! (120061) with the Harps echelle spec- 
trograph. Fourteen oscillation frequencies were detected in the 
power spectrum between 3 and 6 mHz with amplitudes in the 



3.1. Input Physics 

Three different stellar evolution c odes are used for these 
computations: the Gen eva code ((Mevnet & Maedei] l2000t 
Eggenberger et al. 200 7]), the CLES code (Code Liegeois 



d' Evolution Stellaire. IScuflaire et al.ll2007bl) and the C ESAM 
code (Code d'Evolution Stellaire Adaptatif et Modulaire. iMorell 
fl997h . 

These three stellar evolution codes are described in detail in 
the above mentioned references. We s imply note that the Geneva 
code uses the MHD e quation of state dHummer & MihalasS l988: 
lMihalaseta"i]|1988l: iDappenet all [19881) and the NACRE nu- 
clear reaction rates dAngulo et al.lll999|). The CLES cod e uses 
the OPAL equation of state (Rogers & Navfonovl l2002l) in its 
2005 version and the NACRE nuclear reaction rates, while the 
version of the CESAM code used for the present computa- 
tions includes the EFF equation of state (Eggleton e t alj 19731) 
and th e nuclear reaction rates given by Caughlan & Fowler! 
( 1988 j). All codes use th e opacity tables from the OPAL group 
(Iglesias & Rogers 1996) complemented at low temperatures 
with Alexander & Ferguson (1994) opacities and the standard 
mixing-length formalism for convection. The diffusion due to 
the concentration and thermal gradients is also included in the 
computations, but the radiative acceleration is neglected as it is 
negligible for the structure o f low-mass stellar m odels with ex- 
tended convective envelopes dTurcotte et al.ll 1 998b . 

A solar calibration has first been performed with the three 
evolution codes. A mixing-length parameter a = 1.7998 and 
an initial helium mass fraction Y\ = 0.2735 is then determined 
with the physics included in the Geneva code. The input physics 
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used by the CLES code leads to a = 1.8544 and = 0.2755, 
while values of a = 1 .7985 and T; = 0.2848 are found with the 
CESAM code. We note that the higher value of the solar cali- 
brated mixing-length p aramete r found for the CLES code mainly 
results from the use of lKurucz] (1 19981) atmosphere models in the 
CLES code instead of the grey atmosphere approximation in- 
cluded in the other evolution codes. Likewise, the CESAM code 
gives a higher value of the solar calibrated initial helium mass 
faction, which is mainly due to the use of the EFF equation of 
state instead of the OPAL and MHD equations of state intro- 
duced in the CLES and Geneva codes. 



3.2. Calibration methods 

Basically, the calibration of a binary system consists in finding 
the set of stellar modeling parameters that best reproduces all 
observational data available for both stars. The characteristics of 
a stellar model depend on five modeling parameters: the mass 
M of the star, its age (hereafter t), the mixing-length parame- 
ter a = I /Hp for convection and two parameters describing the 
initial chemical composition of the star. For these two parame- 
ters, we choose the initial helium abundance Y\ and the initial 
ratio between the mass fraction of heavy elements and hydrogen 
(Z/X)\. This ratio is directly related to the abundance ratio [Fe/H] 
assuming that log(Z/X) s [Fe/H]+log(Z/X) . In this study, stel- 
lar model s are computed w i th the solar value (Z/X) Q = 0.0245 
given bv lGrevesse" &Noelsld 19931) as well as with the new solar 
mixture of Asplund etal. (2005) complemented with the neon 
abundance of lCunha et alj d2006l) . which results in a lower value 
(Z/X) = 0.0178. 

The binary nature of the system leads to a precise determi- 
nation of the mass of both components. By assuming that both 
components of a binary system share the same age and initial 
chemical composition, three additional constraints are obtained: 
fA = t B , if = Yf and (Z/X)f = (Z/X)f . Consequently, we have 
to determine a set of seven modeling parameters (f, Ma, Mb, 
cxa, as, Yi and (Z/X)i). In the case of 70 Ophiuchi, we also as- 
sume that the value of the mixing-length parameter is identical 
for both components of the system (a^ = <xb)- This assump- 
tion results from the limited number of observational constraints 
available for the B component of the system. Indeed, no astero- 
seismic measurements are available for 70 Oph B and its effec- 
tive temperature is only determined with a large error of 200 K 
(see Table|2]). 

Following the work done for the analysis of the a Centauri 
system, the determination of the set of modeling parameters 
(f, M A , M b , a, T; and (Z/X)i) leading to the best agree- 
ment with the observational constraints is made by using two 
different calibration methods. The first one is based on the 
comp utation of a dense grid of stellar models as explained 
in Eggenberge r et alj d2004bl) . while the second one uses the 
gradient-exp ansion algorithm known as Levenberg-Marquardt 
method (see iMiglio & Monta lban 2005). The reader is refered 
to these two references for a detailed description of the methods. 
In the present study, the grid of stellar models was computed 
by varying the masses of 70 Ophiuchi A and B within their ob- 
servational errors by step of 0.01 M . Concerning the chemical 
composition, the initial helium abundance Y\ was changed be- 
tween 0.240 and 0.290 by step of 0.002. Models of 70 Oph B 
computed with an initial helium mass fraction larger than 0.290 
are indeed found to be too luminous to correctly reproduce the 
observed luminosity. For the initial value of the ratio between 
the mass fraction of heavy elements and hydrogen {ZjX)\, we 
first note that the surface metallicities [Fe/H] s are almost iden- 



tical for models of 70 Ophiuchi with the same initial compo- 
sition and different mixing-length parameters a. Moreover, the 
[Fe/H] s of the models are mainly sensitive to (ZjX)\ and less 
to Y\. As a result, the values of (ZjX)\ are directly constrained 
by the observed surface metallicities; we found that the models 
matching the observed metallicities have (Z/X); ranging from 
about 0.027 to 0.033 (i.e. an initial metallicity [Fe/H]; between 
0.04 and 0.13). The value of the mixing-length parameter was 
changed between 0.8 and 1.3a with a typical step of O.O2a . 
Finally, we note that the age of the models was limited to 15 Gyr 
and that a model of 70 Oph A is considered to have the same 
age as a model of 70 Oph B when the difference is lower than 
0.01 Gyr. 

The theoretical low-^ p-mode frequencies of 70 Oph A 
are calculated using the Aar hus adiabatic pulsation package 
dChristensen-Dalsgaardl IT997h for models computed with the 
Geneva and the CESA M code, while the Liege oscillation code 
dScuflaire et alj l2007al) is used for models computed with the 
CLES code. The value of the mean large separation is deter- 
mined for each model of 70 Oph A by making the average of 
the theoretical separations corresponding to the observed fre- 
quencies (see Sect. 12.4) . A x 2 minimization is then performed 
by defining the following functional that includes the observa- 
tional constraints (classical as well as asteroseismic) currently 
available for both components of 70 Ophiuchi 



X 



= Z 



£<theo _ q> 



-ohs 



(i) 



where N tot is the total number of observational constraints in- 
cluded in the calibration. The vectors C contains therefore all 
observables for both stars: 



C = (M A , M B , L A , L B 



L e A ff ,L e B ff ,[Fe/H]f,[Fe/H]f,Av£) 



The vector C theo contains the theoretical values of these observ- 
ables for the model to be tested, while the values of C obs are 
those listed in Table [2] The vector E obs contains the errors on 
these observations, which are also given in Table|2] 

4. Results 

Using the observational constraints listed in Sect. |2l we perform 
the aforementioned x 1 minimization. In a first time, we fix the 
mixing-length parameter to its solar calibrated value : a a — &b = 
ff . The solar mixture of iGrevesse & Noels d 1993b is used and 
atomic diffusion of helium and heavy elements is included in the 
computation. 

The calibration is first performed by computing a dense grid 
of stellar models with the Geneva stellar evolution code. In this 
way, we find the following solution t = 6.2 + 1.0 Gyr, Ma = 
0.91 + 0.02 M Q , M B = 0.72 + 0.01 M , Y { = 0.266 + 0.015 and 
(Z/X)i = 0.0300+0.0025. The position in the HR diagram of this 
model of 70 Oph A and B (denoted model Mia in the following) 
is given in Fig. [3] The characteristics of this model are reported 
in Table [3] The confidence limits of each modeling parameter 
given in Table[3]are estimated as the maximum/minimum values 
that fit the observational constraints when the other calibration 
parameters are fixed to their medium value. 

In order to compare results obtained by using a different 
stellar evolution code, the same calibration is done by using 
the CLES code and the Levenberg-Marquardt method. This cal- 
ibration (called model Mlb) leads to the following solution: 
t = 6.5 + 1.0 Gyr, M A = 0.90 + 0.02 M , M B = 0.73 + 0.01 M , 
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Table 3. Model s for 70 Oph A and B o btained with three different stellar evolution codes. The Ml models are computed with the 
solar mixture of lGrevesse & Noeli d 19931) and include atomic diffusion of helium and heav y elements. Model M2 is computed with 
the new solar mixture of lAsplund et al.l d2005l) complemented with the neon abundance of lCunha et alj d2006f) and in cludes atomic 
diffusion of helium and heavy elements. Model M3 is computed with the solar mixture of lGrevesse & N oels ( 1993J) but without 
atomic diffusion. These models are computed with a solar-calibrated value of the mixing-length parameter. The upper part of the 
table gives the non-asteroseismic observational constraints used for the calibration. The middle part of the table gives the modeling 
parameters with their confidence limits, while the bottom part presents the global parameters of both stars. 





Model Mia 


Model Mlb 


Model Mlc 


Model M2 


Model M3 




Geneva code 


CLES code 


CESAM code 


Solar mixture: Asp+Cunha 


No atomic diffusion 




70 Oph A 


70 Oph B 


70 Oph A 


70 Oph B 


70 Oph A 


70 Oph B 


70 Oph A 


70 Oph B 


70 Oph A 


70 Oph B 


M/M 


0.89 ± 0.02 


0.73 ±0.01 


0.89 ± 0.02 


0.73 ±0.01 


0.89 ± 0.02 


0.73 ±0.01 


0.89 ± 0.02 


0.73 ±0.01 


0.89 ± 0.02 


0.73 ±0.01 


L/L 


0.53 ± 0.02 


0.15 ±0.02 


0.53 ± 0.02 


0.15 ±0.02 


0.53 ± 0.02 


0.15 ±0.02 


0.53 ± 0.02 


0.15 ±0.02 


0.53 ± 0.02 


0.15 ±0.02 


T cS (K) 


5300 ± 50 


4390 ± 200 


5300 ± 50 


4390 ± 200 


5300 ± 50 


4390 ± 200 


5300 ± 50 


4390 ± 200 


5300 ± 50 


4390 ± 200 


[Fe/H] 


0.04 ± 0.05 


0.04 ± 0.05 


0.04 ± 0.05 


0.04 ± 0.05 


0.04 ± 0.05 


M/M 


0.91 ±0.02 


0.72 ±0.01 


0.90 ± 0.02 


0.73 ±0.01 


0.91 ±0.02 


0.72 ±0.01 


0.90 ± 0.02 


0.73 ± 0.01 


0.90 ± 0.02 


0.73 ±0.01 


r(Gyr) 


6.2 ± 1.0 


6.5 ±1.0 


6.2 ± 1.0 


6.3 ± 1.0 


7.2 ± 1.2 


Yj 


0.266 ± 0.015 


0.266 ±0.018 


0.272 ±0.015 


0.255 ±0.018 


0.254 ±0.018 


(Z/X); 


0.0300 ± 0.0025 


0.0300 ± 0.0025 


0.0300 ± 0.0025 


0.0220 ± 0.0023 


0.0266 ± 0.0025 


L/L 


0.527 


0.152 


0.521 


0.163 


0.535 


0.152 


0.522 


0.164 


0.522 


0.163 


T cS (K) 


5293 


4438 


5283 


4434 


5290 


4395 


5284 


4445 


5286 


4437 


R/Re 


0.865 


0.661 


0.863 


0.685 


0.872 


0.674 


0.864 


0.684 
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Fig. 3. Location in the HR diagram for the models of 70 Oph A 
and B. The dashed lines indicate the boxes delimited by the 
observed luminosities and effective temperatures (with their re- 
spective 1-sigma errors). The dot, the square and the triangle 
correspond to the Mia, Mlb and Mlc models, respectively. The 
open circle indicates the position of the M2 model, while the 
cross corresponds to the M3 model. Note that the symbols cor- 
responding to the Mlb, M2 and M3 models are partly superim- 
posed. In particular, the crosses corresponding to the M3 model 
are superimposed to the squares indicating the location of the 
Mlb model. 



Yi = 0.266 + 0.018 and (Z/X); = 0.0300 + 0.0025. The char- 
acteristics of this solution are given in Table [3] and the loca- 
tion of both components in the HR diagram is shown in Fig. [3] 
A third comparison is finally performed by using the CESAM 
evolution code. The following solution (denoted model Mlc) is 
then found: t = 6.2 + 1.0 Gyr, M A = 0.91 ± 0.02 M G , M B = 
0.72+0.01 M , Yi = 0.272+0.015 and (Z/X); = 0.0300+0.0025. 
The global parameters of 70 Oph AB determined with this so- 
lution are also given in the bottom part of Table [3] while Fig. [3] 
shows the location of both components in the HR diagram. 

From these results, we first conclude that a consistent model 
of the 70 Ophiuchi system that correctly reproduces all obser- 
vational constraints now available for 70 Oph A and B can be 
determined. Indeed, Fig.[3]and Table[3]show that the three solu- 
tions are in very good agreement with all classical observables 
included in the calibration. Moreover, the Ml models correctly 
reproduce the observed value of the mean large separation of 



70 Oph A. Note that the global asteroseismic features of the 
three Ml models of 70 Oph A are very similar; they all ex- 
hibit a mean large separation of about 161.8//Hz with a mean 
small separation between £ — 2 and i — of about 9.5 //Hz. 
As an illustration, the variation in the large and small separation 
6vq2 with frequency for the Mlb model of 70 Oph A is shown 
in Fig. |4] From Table [3] we also conclude that the three differ- 
ent stellar evolution codes give very similar results. We only 
note a slightly higher value of the initial helium mass fraction 
determined with the CESAM evolution code, which is directly 
related to the higher value of this parameter found for the so- 
lar calibration (see Sect. 13. II ). These results are comforting since 
the input physics included in these three codes is similar (al- 
though not striclty identical, see Sect. 13.11 ) and must therefore 
lead to a coherent determination of the modeling parameters of 
the 70 Ophiuchi system. It is also worthwhile to recall here that, 
in addition to the stellar evolution codes, the calibration meth- 
ods were compared. Indeed, the Mia model has been obtained 
by computing a dense grid of stellar models, while the solution 
Mlb has been determined by using an optimization algorithm. 
We thus find that these methods lead to the same determination 
of the global parameters of 70 Oph AB. These results are in good 
agreement with previous computations done for the calibration 
of the a Centaur i system by using the Geneva and the CLES evo- 
lution code (see Eggenberger et al. 2004b; Miglio & Montalban 
l2005h . We also note that the values of the parameters that we 
found for the 70 Ophiuchi system are glo bally in good agr e emen t 
with the results of pr e vious studies by Fernan des et al.l (l998) 
and ICasagrand eet alj ([2007) with smaller error bars resulting 
from the numerous observational constraints now available for 
this system. 

Thus an age of about 6 Gyr with a typical uncertainty of 
1 Gyr is determined for the 70 Ophiuchi system independently 
of the evolution code used for the computation. Concerning the 
initial chemical composition, we find an initial ratio between 
the mass fraction of heavy elements and hydrogen (Z/X); = 
0.0300 ± 0.0025, which is slightly higher than for a solar model. 
This value is directly constrained by the observed surface metal- 
licity [Fe/H] s = 0.04 ± 0.05. Contrary to the initial ratio be- 
tween the mass fraction of heavy elements and hydrogen, a 
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slightly lower value of the initial helium abundance is found for 
70 Oph AB than for a solar model. We however note that the 
error on this parameter is large. Moreover, we recall here that all 
these results have been obtained by fixing the mixing-length pa- 
rameter of both components of the 70 Ophiuchi system to the so- 
lar calibrated value. We are now interested in investigating which 
solutions can be found by relaxing this constraint. 
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Fig. 4. Large and small separations versus frequency for the Mlb 
model of 70 Oph A. The dots indicate the observed values of the 
large separation with an uncertainty on individual frequencies 
estimated to half the time resolution (1.1 //Hz). 



since it mainly depends on the square root of the star's mean 
density. 



-0.2 



-0.3 



O 



-0.4 



-0.5 




-^=0.266 (Model Mia) 
Y,=0.240 



_L 



3.74 



3.73 



3.72 3.71 
Log T efI 



3.7 



Fig. 5. Evolutionary tracks in the HR diagram for two models of 
70 Oph A with different values of the initial helium mass frac- 
tion Y\ and convection parameter a. The model Mia computed 
with Yi = 0.266 has a solar mixing-length parameter, while the 
solution with a lower initial helium mass fraction of T; = 0.240 
exhibits a higher value of this parameter {a = 1.25a ). This 
shows that a decrease in the initial helium mass fraction can be 
compensated by an increase of the mixing-length parameter in 
order to reach the same location in the HR diagram. 



4. 1 . Models with a free mixing-length parameter 

The calibration of the 70 Ophiuchi system is then done by relax- 
ing the constraint of a solar mixing-length parameter. We how- 
ever still assume that the mixing-length parameter of 70 Oph A 
is equal to the one of the B component. We then find that models 
with different values of the initial helium abundance are able to 
correctly reproduce all observational constraints now available 
for the 70 Ophiuchi system. This is due to the fact that a decrease 
(increase) of the initial helium abundance T; can be compensated 
by an increase (decrease) of the mixing-length parameter a to 
match the observed position of 70 Oph A in the HR diagram 
as fo und previously for the calibration of the star Procyon A 
(see [Eggenber ger" et al.l l2005). As an illustration, Fig. [5] shows 
the evolutionary track of 70 Oph A for a model computed with 
an initial helium abundance T; = 0.240 and a mixing-length pa- 
rameter a = \.25a Q . This model shares approximately the same 
location in the HR diagram as the Mia model but exhibits a 
much higher age of about 10.5 Gyr due to its very low initial he- 
lium abundance. We thus obtain a series of models with approx- 
imately the same non-asteroseismic features as those of the Ml 
models computed with a solar calibrated value of the mixing- 
length parameter. Moreover, the mean large separation of these 
models (which is the only asteroseismic quantity included in the 
X 1 functional) is also very close to the value of the Mia model, 



The observed location of 70 Oph B in the HR diagram only 
provides weak constraints to discriminate between these models 
because of the large error on its effective temperature and be- 
cause the model of 70 Oph B spend most of its time within the 
observational box in the HR diagram during its main-sequence 
evolution due to its low mass of 0.73 M . Only solutions with an 
initial helium abundance T; larger than about 0.290 (correspond- 
ing to an age lower than about 2 Gyr) lead to models of 70 Oph B 
with a too high luminosity to reproduce the classical observa- 
tions. This explains why there is a series of models of 70 Oph AB 
with very different ages and initial helium abundances that well 
reproduce the global stellar parameters considered in the cali- 
bration. We however note that these models exhibit different as- 
teroseismic features although they are characterized by the same 
value of the mean large separation. In particular, these models 
have different values of the mean small separation £vo2 since 
their ages differ. This is shown in Fig. [6] where the large and 
small separations of 70 Oph A for the Mia model and for the 
model computed with T; = 0.240 are compared. We see that both 
models are characterized by the same values of the large sepa- 
rations, while the mean small separation <5vo2 of the Mia model 
is significantly larger (3 /iHz) than the one of the model with a 
lower initial helium abundance. This illustrates the importance 
of having a precise observed value of the mean small separation 
of 70 Oph A in order to obtain an independent determination of 
the value of the age, the mixing-length parameter and the initial 
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helium abundance of the 70 Ophiuchi system. This is illustrated 
in more detail in Fig. [7] where the value of the mean small sep- 
arations of 70 Oph A for models of 70 Oph AB computed for 
different initial helium abundances and mixing-length parame- 
ters is shown. As expected, the mean small separation is found to 
significantly decrease with the age. The increase of the age of the 
system is directly correlated to the decrease of the initial helium 
abundance, which is in turn directly related to the increase of the 
value of mixing-length parameter needed to reproduce the ob- 
served location of 70 Oph A in the HR diagram and the observed 
mean large separation. We thus conclude that the observation of 
the mean small separation of 70 Oph A will lift the degeneracy 
between the value of the mixing-length parameter and the initial 
helium abundance by adding a strong observational constraint 
on the age of the system and will therefore enable an indepen- 
dent and precise determination of the mixing-length parameter 
and the initial helium abundance of the 70 Ophiuchi system. 



dance F; between 0.260 and 0.290 when the constraint of a solar 
mixing-length parameter is relaxed. 
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Fig. 6. Large and small separations corresponding to the mod- 
els of 70 Oph A whose location in the HR diagram is shown in 
Fig. [5] The solution with Y\ = 0.266 has a solar mixing-length 
parameter, while the model with Y[ = 0.240 has been computed 
with a = 1.25ar . 

The mean small separation of 70 Oph A has unfortunately 
not been observed. The mono-site nature of current asteroseis- 
mic observations coupled to the low resolution of the time series 
do not enable such a measurement. As discussed in Sect. 12.41 
these observations nevertheless suggest that the value of the 
mean small separation of 70 Oph A is included in a frequency in- 
terval centered around the daily alias of 1 1 .57 //Hz introduced by 
the mono-site observations and delimited by the frequency reso- 
lution of the time-series (2.2 //Hz). As shown in Fig. [7] this con- 
straint on the small separation rules out models of 70 Oph AB 
with an age higher than about 7 Gyr, an initial helium abun- 
dance Yi lower than about 0.260 and a mixing-length parameter 
a larger than about 1.05ar o . We thus conclude that current aster- 
oseismic observations suggest that the age of the 70 Ophiuchi 
system lies between 2 and 7 Gyr with a mixing-length parameter 
between about 0.8 and 1.05or Q and an initial initial helium abun- 
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Fig. 7. Mean small separation i5vo2 of the component A for differ- 
ent models of 70 Ophiuchi. Dots and open symbols correspond 
to solutions with an age lower and higher than 7 Gyr, respec- 
tively. Dotted lines indicate the allowed frequency interval for 
the mean small separation suggested by current asteroseismic 
observations. 



4.2. Models with the new solar mixture 

In order to test the sensitivity of the global parameters deter- 
mined for 70 Oph AB on the adopted solar mixture, we decided 
to redo the whole ca libration using the solar mixture proposed by 
lAsplund et al.l (120051) complemented with the neon abundance 
of ICunha et al. (2006). This results in a lower solar value of the 
ratio between the mass fraction of heavy elements and hydro- 
gen (Z/X) Q = 0.0178. We recall here that solar models evolved 
with this new abundance mixtu re give worse agreement with he- 
lioseismic constraints (see e.g. Guzik 2006). For these compu- 
tations, we use OPAL opacity tables calculated with the corre- 
sponding solar mixture and the mixing-length parameter is fixed 
to its solar calibrated value. We then find the following solution: 
t = 6.3 ± 1.0 Gyr, M A = 0.90 + 0.02 M , M B = 0.73 + 0.01 M Q , 
Yi = 0.255+0.018 and (ZjX\ = 0.0220+0.0023. The position in 
the HR diagram of this model of 70 Oph A and B (denoted model 
M2 in the following) is given in Fig.[3]and its global features are 
reported in Table [3] We note that this solution calibrated with 
the new solar abundances presents a different initial chemical 
composition, but no other significant deviation from the global 
parameters of the Ml mo dels computed with the solar mixture of 
iGrevesse & Noelsl (119931) . Indeed, the Ml and M2 models share 
the same age and the same locations in the HR diagram. A lower 
value of the initial helium abundance is then found for the M2 
model in order to compensate for the decrease in the initial ra- 
tio between the mass fraction of heavy elements and hydrogen 
resulting from the lower solar val ue (Z/Z) . As it happens fo r 
the Sun and for a Centauri A (see Miglio & Montalban 2005), 
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the models of 70 Oph A computed with the solar mixture of 
iGrevesse & Noels! ( 1 1993b exhibit a deeper convective zone com- 
pared with the M2 model, but the uncertainty in the observa- 
tional constraints do not enable to reveal this difference that can 
be masked by other choice of parameters. 



4.3. Models without atomic diffusion 

The effects of atomic diffusion on the calibration of the 70 
Ophiuchi system are also studied by computing models with- 
out diffusion. For these computations, the solar mixture of 
IGrevesse & Noels! fl993) is used and the mixing-length parame- 
ter is fixed to its solar calibrated value. Thus the same calibration 
is made as for the Ml models except for the inclusion of atomic 
diffusion. We then obtain the solution (hereafter refered as M3) 
t = 7.2 + 1 .2 Gyr, M A = 0.90 + 0.02 M G , M B = 0.73 + 0.01 M , 
Yi = 0.254 + 0.018 and (Z/X); = 0.0266 + 0.0025. The char- 
acteristics of this solution are given in Table [3] while the loca- 
tion of both components in the HR diagram is shown in Fig. [3] 
By comparing model M3 with the Ml models, we see that the 
initial chemical composition differs. Indeed, the Ml models are 
characterized by a higher initial ratio between the mass fraction 
of heavy elements and hydrogen than the M3 solution because 
of the inclusion of atomic diffusion. Moreover the initial helium 
abundance of the model without diffusion is lower and this re- 
sults in an increase of the age of about 1 Gyr. 

Figure [8] shows the effects of atomic diffusion on the fre- 
quencies of 70 Oph A. The differences between the frequen- 
cies of the M2 and the Mlb models of 70 Oph A (defined as 
v„o(Mlb)-v„o(M3)) show that the slightly lower value of the sur- 
face abundance in helium of the model computed without atomic 
diffusion (difference of about 4 %) does not significantly change 
the value of the frequencies. For a low-mass star like 70 Oph A, 
the effects of atomic diffusion are indeed small because of the 
large mass contained in its convective zone. We thus conclude 
that the change in the frequencies introduced by the different 
value of the surface helium abundance that implies a different 
opacity and, therefore, a different location of the boundary of 
the convective zone is too small to be clearly revealed. 




v CuHz) 



Fig. 8. Differences of frequencies between the model of 
70 Oph A computed with the inclusion of atomic diffusion (solu- 
tion Mlb) and the model of 70 Oph A computed without atomic 
diffusion (solution M3). 



5. Conclusion 

Our calibration of the 70 Ophiuchi AB shows that a solution cor- 
rectly reproducing all asteroseismic and non-asteroseismic ob- 
servational constraints now available for both stars can be de- 
termined. An age of 6.2 + 1.0 Gyr is found with an initial he- 
lium mass fraction Y\ = 0.266 + 0.015 and an initial metallicity 
(Z/X)i = 0.0300 + 0.0025 when atomic diffusion is included 
and a solar value of the mixing-length parameter is assumed. 
Concerning the mixing-length parameter, we found that current 
asteroseismic measurements do not allow a firm independent de- 
termination of this value that will only be possible with a precise 
measurement of the mean small separation. The effects of atomic 
diffusion and of the choice of the adopted solar mixture on the 
results of the calibration were also studied. We then found that 
a change of the solar mixture results of course in a different ini- 
tial chemical composition but does not change significantly the 
other global parameters. Models computed without atomic dif- 
fusion also have a different initial chemical composition than 
models including diffusion but they also exhibit a slighlty higher 
age (difference of about 1 Gyr). 

We finally mention that we took the opportunity of the cali- 
bration of the 70 Ophiuchi system to compare and test the the- 
oretical tools used for the modeling of stars for which p-modes 
frequencies are detected. Thus, the same analysis was performed 
with three different stellar evolution codes including a coherent 
input physics: the Geneva code, the CLES code and the CESAM 
code. Two different calibration methods were also used: the com- 
putation of a dense grid of stellar models and a minimization 
algorithm. We found that the different evolution codes and cali- 
bration methods lead to perfectly similar and coherent solutions. 
This result is comforting in the perspective of the theoretical 
interpretation of asteroseimic data that are expected from the 
CoRoT and the Kepler space missions and from future ground- 
based campaigns. 
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